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Numerical investigations on nonanalytic quantum kicked systems is presented in this paper. The 
power-law type localization is found to be universal in the nonanalytic systems just like that the ex- 
ponential type in the analytic systems. With increasing the perturbation strength, a transition hom 
perturbative localization, to pseudo-integrable regime, to dynamical localization and to complete 
extension is clearly demonstrated. The diflterent regimes are characterized by the diflterent relations 
between the localization length and the system's parameters. Finally, we investigate the difltusion 
behavior and find that the quantum suppression is relaxed in the nonanalytic systems. 
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INTRODUCTION 

The quantum dynamics of classically chaotic Hamilto- 
nian systems has attracted great interest in recent years. 
Much work has been devoted to study of the periodically 
"kicked" systems. The prime examples are the kicked ro- 
tor (KR)P,|3| and the kicked hydrogen atom In these 
models, an interesting phenomenon is the dynamical lo- 
calization, namely, the quantum suppression of classical 
diffusion. This phenomenon was first discovered numer- 
ically by Casati et al in the KR model , and later on 
confirmed by several experiments such as Rydberg atom 
in microwave field . The deviation from maximal chaos 
in quantum model is entirely related quantum interfer- 
ence effects only. It was shown that the mechanism of 
such a "quantum suppression of classical chaos" is the 
localization of all eigenfunctions in the unperturbed mo- 
mentum space, which similar to the well-known Anderson 
localization in the solid state physics. 

However, these discussions are mainly restricted to the 
analytical systems. On the other hand, nonanalytic sys- 
tems also emerge from the study of more concrete phys- 
ical models, such as, Fermi-acceleration model^], billiad 
system ^ |7| and a kicked particle in ID infinite square 
potential wellQ. Recently, some pioneering works have 
been done towards the systems in which the analyticity 
condition is broken down 0, Hol llll | . Many novel prop- 
erties are found, such as the important role of the can- 
tori structure in the localization. Their discussions are 
mainly towards a special discontinuous model, where the 
derivative of the potential is discontinuous. 

The main purpose of this paper is to observe the quan- 
tum behavior, especially the localization phenomenon in 
the general nonanalytic systems. As will be seen later, 
the nonanalytic systems reveals distinctive properties. 
The power-law type localization is found to be universal, 
just like the exponential localization in analytic models. 
In particular, the power exponent does not relate to the 
derivative order of the potential. It rests on the conver- 
gence rate of the corresponding Fourier series. The lo- 



calization length properly defined in light of the informa- 
tion entropy reveals four different regimes with increas- 
ing the perturbation strength: perturbative localization; 
pseudo-integrable regime; dynamical localization; com- 
plete extension. Each regime is characterized by differ- 
ent relation between the localization length and the sys- 
tem's parameters. Moreover, we investigate the diffusion 
behavior and show that the suppression of the energy 
diffusion can be relaxed in the nonanalytic systems. 

The paper is organized as follows. In Sec. II we in- 
troduce the nonanalytic Hamiltonian considered in the 
paper and discuss its classical diffusion property. The 
power-law type eigenstates is shown in Sec. III. In Sec. IV, 
we calculate the localization length defined in terms of 
the information entropy. The energy diffusion in nonan- 
alytic systems is discussed in Sec.V. Finally, Sec. VI is the 
concluding remarks. 



NONANALYTIC MODELS AND CLASSICAL 
DYNAMICS 

The considered Hamiltonian takes the form (in dimen- 
sionless units). 



where the potential V{q) is 27r-periodic function, taking 
the following form under the Fourier expansion, 

oo 

V{q) = E a„ sin(nq) + 6„ cos(nq). (2) 

n=l 

To simplify our discussions, let fln = 6n = 0, then 



V{q)^V^{q)^C^Yl 



sin(ng) 



(3) 



where the Cq is the normalizing constant so that the 
maximum value of the potential be unit. The parameter 
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a represents the convergence rate of Fourier series. As 
a goes to infinity, only the first term of the series re- 
mains, which is just the sinusoidal function. In this case 
the system is equivalence to the KR model under a dis- 
placement transformation. The potential functions are 
plotted in Fig.l. 
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As well known, the KAM theorem is the footstone of 
modern dynamics theory. It states that invariant tori 
(KAM tori) of irrational winding number will keep up 
under a small perturbation whereas the tori correspond- 
ing to the rational winding number are broken down. 
Fortunately, the broken tori have zero-measure set in the 
phase plane. In the original proof of the theorem, the an- 
alyticity is strictly required. However, it is found later on 
that, the analyticity condition can be relaxed to C'^'^^^'^ 
|l2l | for a twisted map. Applying this result into our sys- 
tem, it means that for a > 6, KAM theorem is available; 
As a < 6, the system is non-KAM system. 




FIG. 1: The nonanalytic potential functions Vaiti) 



Before we go further, let us discuss the symmetry prob- 
lems of the system. Clearly the time-reversal invariance 
is kept for any a. However, for a 7^ cx), two essential 
modifications are made: 1) Vaisi) turns to be nonana- 
lytic (e C^''^), the order of its derivability is /3 = [a] — 2, 
e.g. its (/3 + l)-order derivative diverges or discontinu- 
ous at q = 0. Here [x] represents the minimum integer 
larger than or equals to x; 2) Spatial symmetry is broken 
down, the parity is no longer a good quantum number. 
As shown later, these two changes will substantially alter 
the quantum properties of the system. 

Classically, the Hamiltonian (1) corresponds to a 
twisted map, 

, 914 (q„) 

Pn+l = Vn + K , 

Oqn 

9n+l = (In + TVnW mod{2TT) . (4) 

The classical parameter K — kT determines its classical 
dynamics. In our calculations we fix T = 1. As a be 
an integer, the potential can be expressed analyticly, for 
example. 




FIG. 2: Local phase plane for various parameters. As a = 2, 3 
and k — 0.1, the phase plane is ploted by tracing one trajec- 
tory originated from the stochastic region for ten thousands 
or more periods. 

The main characteristic of the non-KAM system is 
the broken-up of the KAM tori under a small pertur- 
bation so that the diffusion can take place. In Fig. 2, 
we plot the local phase planes of systems corresponding 
to a = 2, 3, oo. It is shown that phase plane become 
more and more chaotic as we increase the perturbation 
strength K. At small value of = 0.1, the phase plane 
of the nonanalytic systems show complicated structures, 
such as the stochastic webs and cantori. The chaotic 
region is connected and the diffusion can occur. As 
K is larger than a critical value Kcr, the phase planes 
become fully chaotic, or ergodic. The critical values 
can be estimated approximatly from the losing stabil- 
ity of the primary resonances. Approximately, it gives 
Kcria = 2) ~ 3.5,Kcr{a = 3) ~ 4.5,Kcria = 5) ~ 4.6 
, 2cn)i Kcr{a = oo) ~ 5.2. In the intermediate perturbation 
K — 2, even though the KAM tori are completely broken 
down, there are some invariant structures, such as small 
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islands, cantori. 

In calculating the diffusion coefficient for a given K, 
we take 1000 points starting from stochastic regions, and 
all the initial trajectories evolve for several thousand pe- 
riods. Averages are taken over 1000 trajectories for each 
time period. It is found that the energy diffusion is asym- 
potically linear for all values of K . The diffusion coeffi- 
cient D{=< AP^ > /n) (n is the time in unit of T) ver- 
sus K is plotted in Fig. 3. As the K is large enough i.e. 
K > Kcr, the diffusion coefficients show D ~ b^K^ in the 
quasilinear approximation, where ba = 0.843, 0.550 and 
0.500 for a be 2,3 and oo, respectively. In this case, the 
phase plane is fully chaotic. Below these threshold values, 
diffusion coefficients show various behavior for different 
a. At a 2, a scaling law like D ^ K^/^ is obviously 
demonstrated. At a — 2>, as K < 0.5, D ^ and 
then a 13 fiT''/^ scaling law. At a = cx), as if > 1 
the diffusion coefficient abruptly rises to the quasi-linear 
approximation regime. 
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X/T is so called quasi-energy, is the eigenstate or 
Floquet state. 

The Floquet states and the quasi-energy can be ob- 
tained by diagonalizing U (T) within a large number of 
plane wave bases |n>,(n = — ^,...^). In our calcula- 
tions N is kept at 1024 and fi = T = 1. The elements of 
matrix are Unm — {"^W {T)\ni) . The fast Fourier trans- 
formation (FFT) is employed to transform the wave func- 
tion between the position representation and momentum 
representation in our numerical calculations. Detailed 
description of the numerical method refer to |^ 0| . 
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FIG. 4: The typical eigenvectors in momentum representa- 
tion for different parameters. The smooth curves represent 
the Lorentz-like distributions (Eq.8). The relation between 
power exponent and the parameter a is 7 = 2a. 



FIG. 3: Classical diffusion coefficients versus the perturba- 
tion strength for various a. 



In the KR model, the eigenstates are always exponen- 
tially localized in the momentum space. For the non- 
analytic systems, the things are quite different. The 
time-reversal invariance imposes certain condition on the 
eigenvector that is. 



(7) 



EIGENSTATES 

Now we turn to the quantum properties of the nonana- 
lytic models and to see how the classical chaos manifests 
itself in the quantum dynamics. The evolution operator 
over the period T of the kick is given by 



U{T) = exp{' 



'^^)exp( )exp{- 



Ah 



)■ (6) 



It is unitary and satisfies following eigenvalue equation, 
C/(T)|$A >= e~Ti"|$A >; Here the eigenphase A is real, 



with a being a constant phase, independent of n. So in 
Fig. 4, only the half number of the components in the 
momentum coordinates are plotted. We find that, for 
a small perturbation (k=0.1), i.e. in the perturbative 
regime, the envelope of Floquet states in momentum rep- 
resentation can be well fitted by Lorentz-like function 
(see left column of Fig. 4): 



Cl 



1 + C2\n — nd''' 
with Cl, C2 are the fitted constants. 



(8) 
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The power-law localization is clearly demonstrated. In 
particular, the power exponent 7 does not relate to the 
P, the derivative order of the potential. It rests on the 
convergence rate of the corresponding Fourier series, i.e. 
a simple relation is given numerically, 

7 = 2a. (9) 

With increasing the perturbation strength, the eigen- 
state shows a chaotic band whose width is proportional 
to k. Outside the band, the tails of eigenvector is found 
to decay, on average, as power-law. The relation between 
the power exponent and the index a expressed by Eq. (9) 
still holds (see right column of Fig. 4) . 

In fact, the power-law type localization of the eigen- 
states can be traced back to the structure of the ma- 
trix U. In the KR model, the values of matrix elements 
Umm+n dccay faster than exponential when n exceeds 
the band width h which is proportional to k, thus the 
elements outside this diagonal band can be safely re- 
garded as zero. Within the band of width 6, the ele- 
ments are proved to be pseudorandom 0. However, in 
our model the situation is different. Careful analysis 
yields that the elements outside the band decay as a 
power law with |t/m,m-i-n| ~ 1/n". We have calculated 
({/^)„(= (C^m,m+ri)) (the average is done over m). The 
typical slope of the curves over a large range is approx- 
imately minus 2a (see Fig. 5). And the diagonal band 
width b in our model is also found to be approximately 
proportional to k (Fig. 5b). This kind of band random 
matrix, describes a new class of physical system, e.g . 
nonanalytic system, has attracted recent attention|l4j. 



LOCALIZATION LENGTH 

As shown above, for nonanalytic system, the quan- 
tum localization is characterized by the power-law decay 
of the eigenvector in momentum coordinates. To give a 
quantitative description of the localization phenomenon 
in nonanalytic svstems, let us investigate the localization 
length advocated in|3. [l5|: 

/ = iVea;p(< H > -Hr), (10) 

Here H denotes the information entropy for a (N-l-1)- 
dimensional eigenvector (^\{n){n — — ■^,...-y), H = 
— $|(n)Zn(<i>^(n)), and the symbol < ... > repre- 
sents an average over all eigenvector; The normalizing 
coefficient Hr determines the maximal value of entropy. 

The above definition is one of commonly used in the 
analytic models such as KR model, it can be extended to 
nonanalytic systems for the information entropy defined 
for a eigenvector is convergence for a > 0.5. 

At a = 00, our model is equivalence to the well studied 
model KR model. In this case, there exist a additional 




FIG. 5: The averaged matrix elements < > for different 
parameters a and k. a) right; b) left. 

integral, the parity. Together with the time-reversal sym- 
metry, the eigenfunction can be expressed by a real func- 
tion with the relation ^\{n) — zL^\{—n). So the total 
number of the independent components of is A^i — N/2, 
it satisfies following Gaussian distribution (e.g. refer to 

Wi\^xin)\) = 2ViVi7^e-*'(")^^ (11) 

Then this analysis gives a implicit expression of the max- 
ium entropy as Hn ~ ln{QA8N). 

In the case a 7^ 00, the parity symmetry is broken, the 
eigenfunctions have different real and imaginary parts. 
Both of them satisfy the Gaussian distribution except 
for A^i — N, then the module of the components should 
satisfy following Wigner-type distribution, 

W{\'fx{n)\) = 27Vi|$;,(n)|e-*'(")^^ (12) 

Our numerical calculations verify the above distribuation 
(Fig. 6) and gives the maximal value of information en- 
tropy as Hr ~ ln{0.7N). 

We calculate the dependence of the localization length 
on the perturbation strength for the different parame- 
ter a (Fig. 7a). It is shown, the plots of the localiza- 
tion length vs. perturbation amplitude k can be divided 
into four parts: perturbative localization regime, pseudo- 
integral regime, dynamical localization regime, and sat- 
uration regime. Each regime is characterized by the dif- 
ferent relation between the localization length and the 
perturbation strength. 

In the perturbative regime (fc < 1), average entropy 
shows linearly proportional to fc, i.e. ln{l) ^ k (Fig. 7c). 
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FIG. 6: The distribution of the typical completely extended 
eigenstates for various a. The smooth curves are distribution 
expressed by Eq.(ll) and Eq.(12), respectively. 

In this regime, the Floquet state in momentum represen- 
tation can be well fitted by the Lorentz-like distribution. 
If k is large enough, one observe a kind of saturation, i.e. 
the localization length almost keep constant, no longer 
increases with the k. This implies that eigenstates have 
become completely delocalized (Fig. 7b). The limit value 
of the length almost tends to N, the maximum localiza- 
tion length. 

In the intermediate k i.e. N > I and fc > 1, the typical 
property of the eigenfunctions is the power-law tail. The 
power exponent tightly relates the convergence rate of 
the Fourier expansion of the external potential. In this 
regime, as the perturbation is large enough, i.e. k > kcr, 
the phase space will show ergodicity. In this case, it is 
clearly shown that the localization length is proportional 
to the square of the perturbation strength i.e. classical 
diffusion rate Dd (Fig.7d, upside of the dashed Hue), 

I - fc2. (13) 

However, in the case of that perturbation is not so strong, 
the phase space seems not ergodic. In this case, even 
though the KAM tori are broken down, the cantori and 
other invariant structures may act as barriers for the 
quantum motions, if the flux through cantori is less than 
one Planck's cell (for example, see middle column of 
the Fig. 2). The eigenstates are expected to localized on 
these invariant structures. Indeed, in such situation the 
quantum system looks as if classical ii rteg rable, i.e. so 
called pseudo-intergrable regime |^llolllll|. In this case, 
the eigenstates still can be approximately fitted by the 



FIG. 7: a) Localization length vs. perturbation strength 

for different a; b) In the saturation regime, a typical com- 
pletely extended eigenvector for a — 5 and k = 300; c) In 
the perturbative regime, it is clearly shown ln{l) ~ k; d) The 
dashed line corresponds to the critical values Kd for differ- 
ent a, approximately. In the downside of the dashed line, i.e. 
pseudo-integral regime, it is clearly shown I ^ k. In the dy- 
namical localization regime (upside of the dashed line), it is 
clearly shown I ^ k^ . 

Lorentz-like distribution (Fig. 8). Furthermore, wc find 
that in this regime the localization length show linearly 
proportional to the perturbation strength (Fig.7d, down- 
side of the dashed line), 

I - k. (14) 

This regime become more and more narrow as we de- 
crease the a. As a < 1.5, the classical coefficient tends 
to infinite which means that flow rather than diffusion 
occurs. In this case, the localization length transits di- 
rectly from perturbative regime to the square scaling-law 
regime (see Fig.7d, case a = 1.2). 

ENERGY DIFFUSION 

The energy diffusion of KR model has been well stud- 
ied. After a typical time scale, the energy diffusion will 
be completely suppressed. For this model, Floquet states 
show the strong localization (exponential decay) in mo- 
mentum representation, which implies that, only a finite 
number of Floquet states can be excited for any initial 
state. The question is that what happens to the energy 
diffusion in the nonanalytic models, where the Floquet 
states are not so 'strong' localized, only the power-law 
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FIG. 8: Typical eigenstates in the pseudo-integral regime. 
The smooth curves represent the Lorentz-like distributions 
(Eq.8). 

type localization is observed from the above sections. In 
the dynamical localization regime, out of a chaotic band, 
the tails of the Floquet states behave as |$a("^)P ^ 
so the average energy approximates to Yin 7^ ^ ^n^, 
which is divergence for a < 1.5. From this analysis, we 
conjecture that the diffusion suppression will be relaxed 
as we decrease the parameter a, i.e. lowering the deriv- 
able degree of the system. 

We make a numerical simulation to test the above de- 
duction. To ensure the precision of the numerical simu- 
lation, both coordinate and momentum spaces should be 
large enough, so a large number (N=32768) of Fourier 
components are used in our computations. We rescale 
the coefficients Ca so that the classical diffusion rate 
Del — Dqi = in the quasilinear approximation. To 
satisfy the condition, = 0.129,0.219,0.759, and 1.00 
at a = 1.52, 1.60,2.00 and oo, respectively. Our results 
are shown in Fig. 9. The complete suppression of the dif- 
fusion is relaxed, instead, a kind of quantum diffusion is 
clearly evident as a be 1.52 and 1.60. Our system is iso- 
lated from the environment, the observed delocalization 
is completely due to the nonanalyticity of the the system, 
which leads to the possibility of a long-range hopping in 
the momentum space. On the other aspect, we find that, 
as a > 2, the diffusion is completely suppressed. 

CONCULSIONS 

In summary, nonanalytic systems show many inter- 
esting properties, the parameter a which represents the 



FIG. 9: The energy diffusion for different parameters a at 
k = a. Classically, the 1000 initial points are distributed in 
p = with the random q and evolved by the map Eq.(4) at 
a = oo. For other a, since the Ca is recaled, in the quasilin- 
ear approximation, the classical evolution shows the similar 
behavior. Quantumly, the initial state is chosen as the ground 
state |0) of unperturbed system. The curves from the up to 
down correspond to a = 1.52, 1.60, 2.00 and oo. It is clearly 
shown that as a > 2 the diffusion is completely suppressed. 
Otherwise, a kind of quantum diffusion occurs. 



convergence rate of the Fourier series determines the 
type of the evolution matrix, therefore plays an impor- 
tant role. Even though our discussion is towards a spe- 
cial choice of the potential, the conclusions are generic. 
The Fourier coefficients of any nonanalytic potential take 
form of c„/n". Previously known nonanalytic models 
are Fermi-acceleration model[3, billiad system 0,01 and 
a kicked particle in ID infinite square potential wellQ. 
Their common feature is the algebraic decay of the evo- 
lution matrix elements, the corresponding parameter a 
equals to one, two and two, respectively. In the gen- 
eral nonanalytic case, our studies reveal four different 
regimes: perturbative, pseudo-integral, dynamical local- 
ization and complete extension. In each regime, the re- 
lation between the localization length and the perturba- 
tion strength is found to follow different scaling law. This 
finding is in agreement with the arguments in previous 
works 9.,KLil] 

On the other hand, it was found that declocalization 
occurs if one let the system interact with an environment. 
This became an active point attracting much attention 
theoretical and experimentally in recent years. Our find- 
ings in this paper suggest that declocalization can be 
realized by breaking the analyticity of the system. We 
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hope our works will stimulate the studies in the direction. 
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Figures Captions 

Fig.l The nonanalytic potential functions Va{q)- 

Fig. 2 Local phase plane for various parameters. As 
a = 2, 3 and k = 0.1, the phase plane is ploted by tracing 
one trajectory originated from the stochastic region for 
ten thousands or more periods. 

Fig. 3 Classical diffusion coefficients versus the pertur- 
bation strength for various a. 

Fig. 4 The typical eigenvectors in momentum represen- 
tation for different parameters. The smooth curves rep- 
resent the Lorentz-like distributions (Eq.8). The relation 
between power exponent and the parameter a is 7 = 2a. 

Fig. 5 The averaged matrix elements < > for differ- 
ent parameters a and k. a) right; b) left. 

Fig. 6 The distribution of the typical completely ex- 
tended eigenstates for various a. The smooth curves are 
distribution expressed by Eq.(ll) and Eq.(12), respec- 
tively. Fig. 7 a) Localization length vs. perturbation 
strength for different a; b) In the saturation regime, a 
typical completely extended eigenvector for a = 5 and 
k = 300; c) In the perturbative regime, it is clearly shown 
ln{l) ^ k; d) The dashed line corresponds to the critical 
values Kci for different a, approximately. In the down- 
side of the dashed line, i.e. pseudo-integral regime, it is 
clearly shown I k. In the dynamical localization regime 
(upside of the dashed line), it is clearly shown I ^ k^. 

Fig. 8 Typical eigenstates in the pseudo- integral regime. 
The smooth curves represent the Lorentz-like distribu- 
tions (Eq.8). 

Fig. 9 The energy diffusion for different parameters a 
at fc = 6. Classically, the 1000 initial points are dis- 
tributed in p = with the random q and evolved by the 
map Eq.(4) at a = 00. For other a, since the Ca is re- 
caled, in the quasilinear approximation, the classical evo- 
lution shows the similar behavior. Quantumly, the initial 
state is chosen as the ground state |0) of unperturbed 
system. The curves from the up to down correspond to 
a = 1.52,1.60,2.00 and 00. It is clearly shown that as 
a > 2 the diffusion is completely suppressed. Otherwise, 
a kind of quantum diffusion occurs. 



